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1 The GUE Conjecture 

1.1 Introduction 

In 1972 H. L. Montgomery announced a remarkable connection between 
the distribution of the zeros of the Riemann zeta- function and the distri- 
bution of eigenvalues of large random Hermitian matrices. Since then a 
number of startling developments have occurred making this connection 
more profound. In particular, random matrix theory has been found to 
be an extremely useful predictive tool in the theory of L-functions. In 
this article we will try to explain these recent developments and indicate 
some directions for future investigations. 

1.2 The Riemann zeta-function 

The Riemann zeta-function is defined by 

oc ^ 

CW = E- {s = a + it,a>l). (1) 



n=l 

It can also be expressed as a product over primes, the Euler product. 



(2) 



for cr > 1. In 1859 Riemann proved that C,{s) extends to a meromorphic 
function on the whole plane with its only singularity being a simple 
pole at s = 1 with residue 1. He further proved that it has a functional 
equation relating the value of ({s) with the value of C(l — s), 

C(s)=xWC(i-s) (3) 

where x(l ~ ■s) = x(s)~^ = 2(27r)~^-r(s) cos(7rs/2). He discovered that 
the distribution of prime numbers is governed by the zeros of C. He was 
led to conjecture that all of the complex zeros p = (3 + i-f of have 
f3 = 1/2. This assertion is the famous Riemann Hypothesis. 

We know that < /3 < 1 for any complex zero of (. Riemann esti- 
mated the zero counting function 

N{T) ^ #{p = /3 + Z7 : < 7 < T} = ^ log ^ + O(logT) . (4) 

Thus, at a height T the average spacing between zeros is asymptotic to 
27r/logT. See [Q for additional background information about 



1.3 Background 



Montgomery |M1] was studying gaps between zeros of the Riemann zeta- 
function in an attempt to prove that the spacings between consecutive 
zeros can sometimes be less than 1/2 of the average spacing. Such a 
conclusion would have led to a good effective lower bound for the class 
number of an imaginary quadratic field. This estimate was not achieved, 
but through the course of his analysis Montgomery was led to conjecture 
that 

1 v-^ ( .sin7ra;^2^ , ^ 

lim — — y 1 = / 1 - ( y dx 5 

where < a < /? are fixed. 

When Montgomery told Freeman Dyson this formula, Dyson re- 
sponded that the integrand was the pair-correlation function for eigenval- 
ues of large random Hermitian matrices, or more specifically the Gaus- 
sian Unitary Ensemble, or GUE. 

The GUE is the limit as ^ oo of the probability space consisting 
oi N X N Hermitian matrices H with a probability measure P{H)dH 
that is invariant under conjugation by any unitary matrix U . Here dH — 
HjXfc d^Hjk Ylj^k d'^Hjk- Mathematical physicists had studied various 
ensembles since the 1950s in connection with work of Wigner in nuclear 



physics. Mehta |Me| has given a thorough treatment of the development 
of the subject. 

Montgomery went on to conjecture that the n-correlation function 
for zeros of ( is the same as that for the GUE; this conjecture came to 
be known as Montgomery's GUE conjecture. 

Odlyzko and Schonhage developed an algorithm that allowed for the 
simultaneous calculation of many values of C(l/2 + it) for t near T in 
average time T*^. This algorithm allowed Odlyzko to do extensive 
computations of the zeros of ^ at a height near zero number 10^''; his 
computations of the pair correlation and nearest neighbor spacing for 
the zeros of C were amazingly close to those for the GUE. His famous 
pictures added much credibility to Montgomery's conjecture. 



1.4 Further evidence for GUE 



In 1996 Rudnick and Sarnak [El£] made some interesting progress on 



the GUE conjecture. To explain their result, number the ordinates of 
the zeros of C{s): < 71 < 72 < . . .. Introduce a scaling 7 = 7^^^ so 
that the 7 have asymptotic mean spacing 1. Then Sarnak and Rudnick 
proved that 



^To^J^ E /fe'---'7jj=/ M/[/,„(x)/(x) dx (6) 



where VFt/_„(x) — WV^„(xi, . . . , x„) is the n-corrclation function for 
the GUE (see section 5.5) and where / is any function satisfying (1) 
/(x + ...,!)) = /(x) for i G M; (2) / is smooth and symmet- 
ric in the variables and decays rapidly as x — > oo in the hyperplane 
Pn '■= {(^^ij • ■ • tXn) ■ J2^j=i — 0}; (3) the Fourier transform /(u) of 
/ is supported in l^il < 2. The condition (1) assures that / is a 

function of the differences of the 7^ . 

Thus, the GUE conjecture has been proven for all correlation func- 
tions, for a limited class of test functions 

1.5 Other L-functions 

The Riemann zeta-function is the prototype for some extraordinary ob- 
jects known as arithmetic L-functions. An arithmetic L-function (or L- 
function for short) has many properties in common with the Riemann 
(■-function. It has a Dirichlet series 

00 

L{s) = ; C^) 

n=l 

it is meromorphic apart from a possible pole at s = 1; it has a functional 
equation 7(s)L(s) — ej{l — s)L(l — s), where e, the sign, satisfies |e| — 
1. L{s) has an Euler product in which the p-th factor is the reciprocal of 
a polynomial in Moreover, an arithmetic L-function is expected 
to satisfy the Riemann Hypothesis, that all complex zeros have real 
part 1/2. (Note that our L-functions are normalized so that the 1/2-line 
is the line of symmetry for the functional equation. See |^ for precise 
definitions.) An L-function is called primitive if it is not the product of 
two L-functions. Primitive L-functions arise from a variety of contexts: 
from Dirichlet characters, from Mellin transforms of certain cusp forms, 
from Galois representations, from algebraic varieties, etc. However, it 
is believed that each primitive arithmetic L-function is associated to a 
cuspidal automorphic representation of GL„i over a number field. 

Rudnick and Sarnak applied their methods to a fairly general L- 
function. They proved the analogue of (6) for any cuspidal automorphic 
L-function over Q, assuming the Riemann Hypothesis for the L-function, 
and with obvious changes to reflect the appropriate scaling of the zeros. 
In particular, the answer did not depend in any way on the distribution 
of the coefficients of the particular L-function. 

The GUE Conjecture is now seen as a universal law governing the 
distribution of zero spacings for all arithmetic L-functions. 

2 Families 



The realization that the compact classical groups of matrices play a role 
in the theory of L-functions arose through seminal work of Katz and 
Sarnak gs|, ^S|. 



2.1 Function Field Analogues 

Katz and Sarnak investigated the distribution of zeros of function field 
zeta-functions. Consider a curve f{x, y) = where / is a polynomial with 
integer coefficients. The zeta-function for / over a finite field can be 
obtained in a simple way from the generating function of the numbers Nn 
of points on that curve in the finite field extensions F^n , n = 1,2,.... It 
is known that this zeta-function is a rational function whose numerator 
is a polynomial with integer coefficients, degree 2g where g is the genus 
of the curve, and that it satisfies the Riemann Hypothesis that all zeros 
have modulus l/y/q- One can order these zeros on the circle \z\ = in 
terms of their angles measured from the positive axis. One then considers 
statistics of the angles. 

Katz and Sarnak proved that, after proper normalization, the n- 
correlation function of these angles, as g and q tend to infinity, is exactly 
the ri-correlation of GUE. They also proved that the nearest neighbor, 
or consecutive spacing, statistic for these angles is the same as that for 
GUE. 

The method of proof involved working with subgroups of \J{N), the 
group of N X N unitary matrices with the Haar measure (see sections 

5.1 and 5.2). It had been shown by Dyson that the n-correlation and the 
spacing statistics are the same for GUE as for the appropriately scaled 
limit as ^ oo for U{N) with Haar measure. (In the physics literature 
the limit of U{N) is called the Circular Unitary Ensemble or CUE; see 

The function field zeta-functions that Katz and Sarnak were work- 
ing with were characteristic polynomials of matrices from subgroups of 
U(A^). To deduce their result they applied a theorem of Deligne about the 
equidistribution of function field zeta-functions among the characteristic 
polynomials of conjugacy classes in these subgroups. 

2.2 Introduction of Families 

Katz and Sarnak investigated the robustness of their theorem. They con- 
jectured that the conclusion remains true under the weaker hypothesis 
that q is held fixed and g oo. On the other hand, they provided ex- 
amples of sequences of curves of increasing genus for which the spacing 
statistic is not the GUE consecutive spacing statistic. 

They also considered special families such as (a) curves of the form 
2/^ = f{x) with all squarefree, monic /, and (b) quadratic twists of 
curves of the form = x{x — l){x — t); that is Ay"^ = x{x — l){x — t) 
where Z\ is a squarefree, monic polynomial over Fg. When q and the 
degree of Z\ — > oo they discovered that the n-correlation and spacing 
statistics again matched up with CUE. However, the reason in each case 
is different depending on the geometric monodromy group. In case (a) 
they computed that the monodromy group was the symplectic group 
and in case (b) the orthogonal group. 



We denote by USp(Af ) (if N is even) the unitary symplectic matrices, 
and by O(A^) the orthogonal matrices. These are subgroups of U(A^) and 
are equipped with their own Haar measure. 

Katz and Sarnak proved that the n-correlation statistic and the spac- 
ing statistic for the hmits of U(A^), USp(A^), and O(A^) are ah the same 
as that for GUE. 

By contrast, the Circular Orthogonal Ensemble (COE) and the Cir- 
cular Symplectic Ensemble (CSE), which are well-known ensembles in 
mathematical physics (see |Me|) have the same underlying symmetry 
groups as O and Sp, but have a measure different from the Haar mea- 
sure. They have different n-correlation and spacing statistics than those 
of GUE and CUE. 

From now on we will use the letters O, Sp, U when referring to a 
statistic associated with the limits of O(A^), USp(A^), and U(A^). These 
statistics are computed for N x N matrices of the appropriate subgroup 
of U(7V) with their respective Haar measure and then N ^ oo with an 
appropriate scaling limit in order to determine the statistics (see 5.4, 5.5, 
and 5.6). 

Since orthogonal and symplectic matrices have eigenvalues which oc- 
cur in complex conjugate pairs, it is clear that the eigenvalue 1 plays a 
special role for these matrices, whereas it does not for unitary matrices. 
Indeed, there are statistics which differentiate O, Sp, and U. In partic- 
ular, the eigenvalue nearest to 1 is such an example. It turns out that 
this statistic is different for all three symmetry types. More generally the 
j-th eigenvalue nearest to 1 is a statistic that is dependent on the group 
(for a summary of the four statistics of interest to us, see section 5.3). 

Another statistic that is different for all three symmetry types is 
"level-density" . The n-level density function is obtained from summing 
a test function at n-distinct eigenvalues. For comparison purposes, note 
that the n-correlation function is obtained by summing a test function of 
n-variables which is a function depending only on the differences of the 
arguments at n distinct eigenvalues. See sections 5.5 and 5.6 for more on 
these statistics. 

Katz and Sarnak showed that for their special families of function 
field zeta-functions these new statistics (j-th nearest eigenvalue to 1 and 
n-level density) match with the appropriate statistics from Sp and O, 
which they had computed. 

In general, they found that if they could compute the geometric mon- 
odromy group associated with a family of function field zeta-functions, 
and if that monodromy group was U, Sp or O, then all four of the statis- 
tics we have been discussing for the function field zeta-functions could 
be proven to match with the appropriate statistic from U, Sp or O. 

Thus, we say that the symmetry type for the family of curves of the 
form = f{x) the symmetry type is Sp; and for the family of quadratic 
twists of — x{x — l){x — t) the symmetry type is O. 

Katz and Sarnak also have examples of families of function field zeta- 
functions where the symmetry type is U. 



2.3 L-Functions Over Number Fields 



Katz and Sarnak speculated that their results for function field zeta- 
functions would have implications for L-functions over number fields. 
Two collections of L-functions which present themselves as natural ana- 
logues to the families above are (a') the collection of all Dirichlet L- 
functions L(s, Xd) where Xd{n) is a real primitive Dirichlet character with 
conductor \d\ (so that d runs through the set of fundamental discrimi- 
nants of quadratic number fields) and (b') the collection of L-functions 
associated with primitive Hecke newforms of a fixed weight. The analogy 
with the function field zeta-functions is as follows. The zeta-function for 
a member of the Sp family (a) above is obtained by counting solutions 
to the equation — f{x) over a finite field. These solutions are counted 
in a field with p elements by the sum X)a=i Xp(/(o)) where the real 
Dirichlet character Xp modulo p appears. This suggests that the family 
(a') could have a symmetry type Sp. Similarly, the zeta-function of a 
member of the O family (b) above is the reciprocal of the p-th factor in 
the Euler product for the L-function of the elliptic curve defined through 
twists of the equation = x{x—l){x — t), the Legendre family of elliptic 
curves. These L-functions are known (by the solution of the Taniyama 
- Shimura conjecture) to be associated to primitive Hecke newforms of 
weight 2. Thus, the collection (b') could well be a family with symmetry 
type O. 

The first evidence that (a') is a family with symmetry type Sp came 
from Michael Rubinstein's thesis |Q . He computed the lowest lying zero 
of L(s, Xd) and the data matched well with the eigenvalue nearest 1 
for symplectic matrices. He examined theoretically the n-level density 
of zeros and showed (for test functions with restrictions on the support 
of their Fourier transforms) that the n-level density functions for the 
zeros of L(s, Xd) are identical with the n-level density functions for Sp. 
Precisely, we index the ordinates of the zeros of L{s,Xd) as 



0<li''^<li'^<... (8) 



and scale using 



ZTT 

Let D* = X]|(i|<r> ^- Then Rubinstein proved that 

E E /(7jf\---,7S')- / W^5p,„(x)/(x)dx, (10) 

as D ^ (XI where Wsp,n is as in section 5.5, provided that / is a Schwarz 
function such that the support of f{u) is contained in 1%! < 1- 

Wsp,n is called the n-level density function for the symplectic group. 

Rubinstein also found evidence that another family has a symme- 
try type O. To describe this family let A{z) — T{n)e{nz) where 



T is Ramanujan's tau-function. It is well-known that Zi is a primitive 
Hecke newform of weight 12 and level 1. The family (c') is obtained 
from quadratic twists of the L-function associated with A namely 



If d < then this L-function automatically vanishes at s = 1 /2 because 
the associated functional equation occurs with e = — 1. Rubinstein com- 
puted the lowest zero of L{A, s, Xd) for > and the lowest zero above 
the real axis for L{A, s, Xd) with li < 0. At this point, wc should men- 
tion that an orthogonal matrix has determinant or —1. So there are 
actually two symmetry types 0+ and 0~ (and the statistics of O are an 
average of the statistics of these). Katz and Sarnak had computed the 
statistics (neighbor spacing, correlations, density, eigenvalue nearest 1) 
for these two symmetry types as well. Rubinstein found that the lowest 
lying zero of L-functions in (c') with d > followed 0+ while the lowest 
lying zero above the real axis for L-functions with d <0 followed 0~. 

Rubinstein considered more generally the twisting of an arbitrary 
arithmetic L-function L{f, s) associated with an automorphic form / on 
GLm by quadratic characters Xd- The results here divide into three cases 
which have to do with signs of the functional equations. If the L-function 
L(/, s,sym^) associated with the symmetric square of / is entire, the 
functional equations for L[f,s,Xd) will always have e = and the 
average over d is exactly as above with the n-level density function being 
Wsp.n - In the case that the symmetric square L{f, s, sym^) has a pole (at 
s = 1) then the sign of the functional equation for L{f,s,Xd) is +1 for 
even characters Xd and —1 for odd characters Xd- Rubinstein averages 
over these two cases separately and discovers that the first case yields an 
n-level density function Wo+,n and in the second case an n-level density 
function Wo-,n (see section 5.5). In each of these cases a restriction is 
placed on the support of the Fourier transform of /. 

2.4 The Diagonal Terms 

Rubinstein's work gave impressive confirmation of the theory, but the 
severe restriction on the support of the Fourier transform is worth investi- 
gating. In fact, this restriction occurs right at the place where something 
interesting is happening with the Fourier transfom of the density func- 
tion. Returning for a moment to the case of the Riemann zeta-function, 
we can illustrate this idea. 

Montgomery's original theorem involved 




(11) 



F{a,T) = 



1 



T*"(t-t')w(7 - V) . 



(12) 



N{T) 



0<7,7'<T 



Here w is a weight function that concentrates at the origin. (Montgomery 
used w{u) = 4/(4 -|- w^).) Assuming the Riemann Hypothesis he showed 



that 

F(a,T) =r-2"(i+o(l)) + |a|+o(l) (13) 

uniformly for — l + e<Q;< 1 — e for any e > 0. 
By the definition of F, 



1 f°° 
mTt) ^ t{i - ~ l') ^ 1 f{a)F{a,T) da. 



(14) 



Thus, Montgomery's theorem gives information about the average be- 
havior of differences between the zeros for test functions r with the sup- 
port of f contained in (—1, 1). Montgomery went on to conjecture (based 
on considerations of the behavior of prime pairs) that for |a| > 1 one has 
F{a,T) = 1 + o(l); this assertion implies (1). Thus, F is not differen- 
tiable at a = 1. In the proof of Montgomery's theorem (via the explicit 
formula and the mean value theorem (23) for Dirichlet polynomials) for 
\a\ < 1, the main term of F arises from the "diagonal" contributions 
of the mean square of a Dirichlet polynomial (i.e. the terms m = n in 
the integral amJhii'm / nY* dt). For a > 1, the off-diagonal terms (i.e. 
m n) contribute to the main-term. 

A similar situation arises in the work of Rudnick and Sarnak and in 
the work of Rubinstein. All of the proofs of these theorems are valid only 
in the range where the diagonal terms dominate. 



2.5 Beyond the Diagonal 

So far we have seen that the theory of families is confirmed by numerical 
data as well as theoretical data up to the diagonal. Bogolmony and 
Keating gave a heuristic derivation of all of the GUE conjecture (i.e. 
all the n-level correlations) based on Hardy-Littlewood type conjectures 
for pairs of primes and pairs of almost primes; this work shows how the 
off-diagonal terms potentially contribute. 



Ozliik |0z| proved an analogue for all primitive Dirichlet characters 



for the pair correlation theorem (12); he obtained a result for |a| < 2. 



This was the first example of going beyond the diagonal. See also [0S| 
and ^ 



It is of great interest that Iwaniec, Luo, and Sarnak ILS | have suc- 
ceeded in going beyond the diagonal in several examples which represent 
three of the symmetry types. They work with the 1-level density func- 
tions assuming only that the Riemann Hypothesis holds for all of the 
L-functions in question. The 1-level density functions may be obtained 
from Wn by taking n = I. They are 

W{0){x)^l + ^So{x), (15) 



(16) 



(17) 



mSp)(x) = l-^^, (18) 

W{\J){x) = l. (19) 

The families considered by Iwaniec, Luo, and Sarnak are related to mod- 
ular forms. Let H^{N) denote the set of holomorphic newforms / of 
weight k and level N. Let (N) denote the weight k level A'' newforms 
for which the associated L-function has a + in its functional equation, 
and H^{N) is the subset of / for which L{f,s) has a — in its func- 
tional equation. Let M*{K, N) be the union of the H*{k, N) ioT k < K 
and similarly define M+ and M~ . They consider the low lying zeros of 
L[f, s) as / varies through one of the M-sets. The average spacing for 
all the zeros of all the L{f,s) with / e H*{k,N) up to a fixed height 
to is asymptotic to 27r/ log(fc^A^). Let be a test function which is even 
and rapidly decaying. They proved that if the support of cj) is contained 
in (-2,2), then 

1 ^ fjflogPN 



\M*(K,N)\ 

-If 



E 4>('-^^)^ rmwm^)dx. (20) 



as KN — > 00. Similar statements hold with M* replaced by M+ and 
M~ and O replaced by 0+ and 0~. 

It should be pointed out that the Fourier transforms of the density 
functions W{0){x), W{0+){x), and W{0'){x) all agree in the diagonal 
range; so it is only when one goes beyond the diagonal that the distin- 
guishing features of these three symmetry types becomes apparent. 

Iwaniec, Luo, and Sarnak also consider the symmetric square L- 
functions of the / e M* and verify that the above statements hold 
with symmetry type Sp and the support of (p in (-3/2,3/2). Also, the 
average zero spacing is 27r/ log(fc^A^^) so N should be replaced by N"^ in 
the argument of (f) in the left hand side of (20). 

Thus, the theoretical and numerical evidence that the zeros of families 
of L-functions depend on the symmetry type of the family is pretty 
convincing. 



3 Moments 



So far we have seen that the eigenvalues of matrices from unitary groups 
are excellent models for zeros of families of L-functions. Now we want to 
take the matrix models a significant step further and argue that the char- 
acteristic polynomials of these matrices on average reveal very important 
features of the value distribution of the L-functions in the family. 



3.1 Moments of the Riemann zeta-function 



We first look at the situation of the Riemann zeta-function and its mo- 
ments. 

To give some background, we cite the theorem of Hardy and Little- 
wood: 



and the theorem of Ingham: 



1 

TJo 



27r2 



(21) 



(22) 



The asymptotics of no other moments (apart from the trivial 0-th mo- 
ment) are known. In general it has been conjectured that 



(23) 



The basic tools for investigating mean-values in i-aspect are the 
mean value theorem for Dirichlet polynomials (due to Montgomery and 
Vaughan) : 



•^0 n<N 



dt= J2 (T + Oin))\an\ 



(24) 



n<Af 



and some sort of formula expressing the function in question in terms 
of Dirichlet polynomials (such as an approximate functional equation) 
such as 

c(.)'^=E^+x(.)'=E^+^(^)' (25) 



n<T^ 



where E{s) should be small on average, s ~ 1/2 + it, t = y/t/{2Tr), x(s) 
is the factor from the functional equation (3), and where 



dk{n) 



(a>l) 



(26) 



so that C('S)'^ is the generating function for dk{n). Note that the mean- 
value theorem for Dirichlet polynomials detects only diagonal contribu- 
tions. 

Conrey and Ghosh |CG2| gave the moment conjecture a more precise 
form, namely that there should be a factorization 



Cfc 



r(l + fc2) 



where 



oo 



n s 



(27) 



(28) 



is an arithmetic factor and gk, a geometric factor, should be an integer. 
Note that by the mean-value theorem for Dirichlct polynomials it is not 
difficult to show that 



f 

T 



E 



dk (n) 



dt ■ 



r{i + k'^) 



(29) 



provided that x = o{T). Thus, an interpretation of gk is 

J^\CHk + ^t)\' dt 



9k 



lim 



1<T „l/2 + it 



(30) 



dt 



assuming that the limit exists, so that gk represents the 'number' of 
Dirichlet polynomial approximations to Ci^)'' of length T required to 
measure the mean square of ({s)''. 

In this notation, the result of Hardy and Littlewood is that gi — 1 
and Ingham's result is that g2 = 2. 

These results can be obtained essentially from the mean- value theo- 
rem for Dirichlet polynomials. To go beyond the fourth moment requires 
taking into account off-diagonal contributions. Goldston and Gonek |GG| 
describe a precise way to transform information about coefficient corre- 
lations X]n<x a('^)i(?^ + r) into a formula for the mean square of a long 
Dirichlet polynomial X]ri<a: a(^i)"'~* where x is bigger than the length of 
integration. 



Using Dirichlet polynomial techniques Conre y and Ghosh | CGI | con- 
jectured that .93 = 42 and Conrey and Gonek | CGo| conjectured that 
(74 = 24024. Meanwhile, Keating and Snaith |KeSnl| computed the mo- 
ments of characteristic polynomials of matrices in U{N) and found that 
for any real x and any complex number s, 



Mu,n{s) 



N 



U{N) 



I det{A ~ Iexp{-ix))f' '^^ = 11 



r{j)r{j + 2s) 



(31) 

where dA denotes the Haar measure for the group U (N) oi NxN unitary 
matrices. To do this calculation, they made use of Weyl's formula for the 
Haar measure (see section 5.3) and Selberg's integral (see section 5.6). 
They also showed that 



lim 



Mjv(s) G(l + s)2 



G{l + 2s) ' 



(32) 



where G(s) is Barnes' double Gamma-function which satisfies G(l) = 1 
and G{z + 1) = r{z)G{z). Note that for s = fc an integer, 



G(l + fc)2 



fe-i 

n 



G(l + 2fc) fj-^{j+n)\ 



(33) 



For fc = 1,2,3 the above is l/l!,2/4!, 42/9! in agreement with the 
theorems of Hardy and Littlewood, and Ingham and the conjecture of 
Conrey and Ghosh. Keating and Snaith argued that one should thus 
model the moments of the zeta-function from to T by moments of 
characteristic polynomials of unitary matrices of size N ~ logT. (More 
precisely, one should take N to be the integer nearest to log^). They 
then conjectured that 

^^-'^'^■lij-^, (34) 

for integer k. The initial public announcements of the conjectures of 
Conrey and Gonek (that — 24024) and of Keating and Snaith ( gk for 
all real k > —1/2) occurred at the Vienna conference on the Riemann 
Hypothesis only moments after it was checked that the Keating and 
Snaith conjecture does indeed predict that g^ — 24024. 



3.2 Moments of L-functions at 1/2 



Subsequently, Conrey and Farmer |CF[ | analyzed known results for mo- 
ments of L-functions at 1/2 and made a general conjecture. (The se m o- 



ments had been considered by a number of authors; see especially |GV|.) 
The conjecture has the shape 

— y L{f, 1/2)''^ , ■'^'"^ , riogX)^^^ (35) 

for some a^, g^, and q{k) where J- is a family of / parametrized by 
the conductor c(/), and X* = J2c{f)<x ^- "^^"^ observations of Conrey 
and Farmer were that gk and q{k) depend only on the symmetry type 
of the family, and that ak depends on the family itself, but is explicitly 
computable in any specific case. Thus, the conjecture is that q{k) = k^ 
for a unitary family, q{k) = k{k+l)/2 for a symplectic family, and q{k) — 
fc(fc— 1)/2 for an orthogonal family. The values of gk were left unspecified, 
but were then predicted as before from random matrix theory by Keating 



and Snaith [KeSn2| and independently by Brezin and Hikami pH|| by 
computing moments of characteristic polynomials of matrices from 0{N) 
and from USp(2iV). Each is a quotient of products of Gamma- functions. 

Thus, for the family of Dirichlet L-functions L{s, Xd) with a real 
primitive Dirichlet character Xd modulo d we have the following results 
(D* = T,\d\<D 1): Jutila proved that 

^Y.^(i^^d}-aMDi) (36) 

\d\<D 



and 



\d\<D 



Soundararajan |3ol] showed that 
1 



D* 



|d|<_D 



6! 



and conjectured that 



i\<D 



04 

10! 



(38) 



(39) 



The general conjecture coming from random matrix theory (see sec- 
tion 5.7), which agrees with the above, is: 



(40) 



\d\<D 



where 



e=i 



(41) 



An example of an orthogonal family where several moments are 
known arises from J^q the set of primitive cusp forms of weight 2 and 
level q {q prime). Then, from results of D uke Duke, Friedlander, 
and Iwaniec | DFI | , Iw anie c and Sa rnak ISl , and Kowalski, Michel, and 
VanderKam |1KMV1[| and |KMV2| , we have 



J2 i(l/2,/)~ai 



J2 L\\/2J)^2a2\ogq^ 



3 „i 



6 „i 



where oi — C(2), 



a4 



a2 = C(2)'n(l + ^ 

p 

a3 = C(2)^n(l-i) (l 

p 

c(2)^n(i-i)'(i+f+ 



(42) 
(43) 
(44) 
(45) 

(46) 
(47) 
(48) 



We have not found a simple expression for ak, though it can be deter- 
mined expUcitly for each k. As before, a general conjecture is: 



1 



fc-i 



£\ 



^ L'=(l/2,/)^2'=-i[] 



(49) 




We believe that in formulas for moments of L-functions over a family 
the power of the log of the conductor and the value of should only 
depend on the symmetry type of the family and that the value of Ofe will 
depend on the family but can always be determined explicitly. 

4 Further directions 

In this section we mention some questions where further research is de- 
sirable. 

4.1 Full moment conjecture 

What are the lower order terms in the moment formulae for |C(l/2+ii) 
and for L(l/2)'^? These are known in a few instances (see Q, |cj for the 
second and fourth moments of C(s)) but not in general. The difficulty is 
that random matrix theory does not "see" the contribution of the arith- 
metic factor Ofe . Lower order terms will likely involve a mix of derivatives 
of afc and secondary terms from the moments of the characteristic poly- 
nomials of matrices. In general, a better understanding of how ^(s) is 
modeled by a characteristic polynomial of a certain type of matrix is 
needed; how do the primes come into play? Perhaps we should think of 
C(l/2-|-it) as a partial Hadamard product over zeros multiplied by a par- 
tial Euler product. Perhaps these two parts behave independently, and 
the Hadamard product part can be modeled by random matrix theory. 

4.2 Distribution of Values 

Keating and Snaith compute explicit formulas for the s-th moment of 
the characteristic polynomials of matrices from O(A^), Sp(A^), and U(A^). 
Consequently the value distributions for these characteristic polynomi- 
als can be explicitly computed; they involve the Fourier transform of the 
s-th moment. Preliminary investigations indicate that there is a good fit 
between the random matrix formulae and numerical data. One particu- 
larly interesting feature of this investigation involves the understanding 
of zeros which occur exactly at 1/2. These seem to occur only for L- 
functions in an orthogonal family. If the sign of the functional equation 
is —1, there is automatically a zero at 1/2. The interesting situation is 
when the sign is -I- and there is still a zero; for example if E is an elliptic 
curve defined over Q then examination of the distribution of values of 
0+ suggest that twists L{E^ s, Xp) of the L-function by quadratic char- 
acters seem to vanish for about X'^/'^{logX)~^^^ values of \p\ < X with 
sign -|-1 in the functional equation. 



4.3 Extreme Values 



How large is the maximum value of |C(l/2 + for T < t < 2T1 It 
is known that the Riemann Hypothesis implies that the maximum is at 
most exp(clogT/ loglogT) for some c > 0. It is also known that the 
maximum gets as big as exp(ci(logr/ loglogT)^/-^) for a sequence of 
T ^ oo for some ci > 0. It has been conjectured that the smaller bound 
(the one that is known to occur) is closer to the truth. However, the new 
conjectures about moments suggest that it may be the larger. 

This question has a number of equivalent and analogous (for L- values) 
formulations. How big can Sir) := i argC(l/2 + iT) be? Assuming the 
Riemann Hypothesis, it is known that S{T) <^ log T/ loglogT but that 
infinitely often it is bigger than c(logT/loglogT)^/^ for some c > 
Which is closer to the truth? What is the maximum size of the class 
number of an imaginary quadratic field (as a function of the discrimi- 
nant) (see |Sh]for a discussion and numerical investigation.)? How big can 
the least quadratic non-residue of a given prime p be (logp or (logp)^"*^? 
See |VI2|] for a discussion of this question? What is the maximal order of 
vanishing of an L- function at 1/2? In terms of the conductor A^, can it 
be as big as logiV/ log log or is it at most the square root of that, or 
something entirely different? All of these questions are related, at least 
by analogy, and they may all have similar answers. It would be interest- 
ing and surprising if in each case it is the larger bound which is closer 
to the truth. 



4.4 Zeros of C'(s) 

Can one use random matrix theory to predict the horizontal distribution 
of the real parts of the zeros of C'? It is known that the Riemann Hy- 
pothesis is equivalent to the assertion that each non-real zero of C'('S) has 
real part greater than or equal to 1/2. Moreover, if such a zero has real 
part 1/2, then it is also a zero of C(s) (and so a multiple zero of C(s))- 
These assertions are the point of departure for Levinson's work on zeros 
of the Riemann zeta-function on the critical line. It would be interesting 
to know the horizontal distribution of these zeros; in particular what 
proportion of them with ordinates between T and 2T are within a/logT 
of the 1/2-line? 

In a similar vein, the Riemann ^-function is real on the 1/2-line and 
has all of its zeros there (assuming the Riemann Hypothesis). It is an 
entire function of order 1; because of its functional equation, ^(l/2-|-iy^) 
is an entire function of order 1/2. It follows that the Riemann Hypothesis 
implies that all zeros of ^'(s) are on the 1/2-line. Assuming this to be 
true, one can ask about the vertical distribution of zeros of f'(s), and 
more generally of ^^"^(5). It seems that the zeros of higher derivatives 
will become more and more regularly spaced; can these distributions be 
expressed in a simple way using random matrix theory? 



4.5 Long Mollifiers, Local Integrals, and GUE 



David Farmer |F1|, |F2| has made two very interesting conjectures having 
to do with ({s). The first is a conjecture about the mean square of ({1/2+ 
it) times an arbitrarily long moUifier. A mollifier is a Dirichlet polynomial 
with coefficients equal to the Mobius /i-function times a smooth function. 
The length of the mollifier is the length of the Dirichlet polynomial. He 
has also conjectured that 

1 r + u)C{l -s + v) iu-a)iv-b) 

TJo C(3 + a)C(l-.s + 6)^* {u + v)ia + b)^' ^'^""^ 

where a,b,u,v are complex numbers with positive real part, and s — 
1/2 + it. These two conjectures are essentially equivalent and imply cer- 
tain parts of the GUE conjecture. It would be interesting to generalize 
these and relate them to the full GUE conjecture. 



5 Appendices 

5.1 The Classical Groups 

— The unitary group U(iV) is the group of N x N matrices U with 
entries in C for which UU* = I where U* denotes the conjugate 
transpose of U, i.e. if U = (uij), then U* = iuJJ)- 

— The orthogonal group 0{N) is the subgroup of U(iV) consisting 
of matrices with real entries. 

— The special orthogonal group SO(Af). This is the subgroup of 
0(iV) consisting of matrices with determinant 1. S0(2iV) leads to 
the symmetry type we have called O"*" and S0(2iV+l) leads to the 
symmetry type we call 0~. 

— The symplectic group USp(27V) is the subgroup of U(27V) of ma- 
trices U for which U ZU^ = Z where denotes the transpose of U 
and 



5.2 The Weyl Integration Formula 

The N X N unitary matrices can be parametrized by their N eigenval- 
ues on the unit circle. Any configuration of N points on the unit circle 
corresponds to a conjugacy class of U(A^). If f{A) = f{6i, . . . On) is a 
symmetric function of N variables, then Weyl's formula gives 




where dA is the Haar measure. Similarly, on Sp(2iV) and S0(2A^) we 
have respectively 

2jy2 n n 

dA = — Y[ (cos TT9j - cos TT9k f Yi ^^^^ 11 ^ (^^) 

j<k j=l 3 = 1 

dA = J\{C0S TT9j - cos TT0kf Yl '^^3 ' (^4) 

i<3 j=l 



5.3 Four Statistics 

Suppose we have a sequence T of -tuples of numbers Tjv = {ti, t2, • • • , ^tv} 
where t\ < t2 < ■ ■ ■ < In such that for each set the average spacing 
tjj^i — tj is asymptotically 1. We write U^n in place of U if we need to 
indicate that U € Tjy. 

— The n-level density of T is W{x) = W{xi, . . .a;„) means that 

lim J2 fitn^---tiJ= [ /(x)W^(x)dx. (55) 

(il,...i„),i^-<N "'K 

for a Schwarz-class /. The sum is over n-tuples with distinct entries. 

— The j-th lowest zero density is vj {x) means that for a test func- 
tion / 

1 f°° 
N^ooN ^ " / f^^^^i^^^ 

n<N 

— The consecutive spacing density is iJ,{x) means that for a test 
function f{x) we have 

1 f°° 

JV— >oo iV — ' In 

i<N-l •'^ 

— The n-correlation density is V{xi, . . . ,Xn) means that for test 
functions / that are symmetric in all of the variables, depend only 
on the differences of the variables (i.e. f{xi + u, . . . , Xn + u) = 
f{xi, . . . , Xn) for all u), and are rapidly decaying on the hyperplane 
Pn : {(-Ti, . . . ,Xn ■ JZ^i — havc, as N ^ oo, 

1 f{tl,---tn)^ [ fix)V{x) dXi... dXn-1 (58) 

as N oo. The spacing and n-correlation densities are universal, i.e. 
the same for each of O, Sp, and U, whereas the n-level and j-th lowest 
zero densities depend on the symmetry type. 



5.4 Gaudin's Lemma 



Associated to each N x N unitary matrix A arc its N eigenvalues e{9j) 
where < < . . . < < 1. We integrate a function F{A) over U(A'') 
by parametrizing the group by the 9i and using Weyl's formula to convert 
the integral into an A^-fold integral over the Oi. 

Often one wants to integrate with respect to Haar measure over U(A^) 
a function F{A) = f{A) ^ f{9i, ...,0n) ofN variables that is "lifted " 
from a function f oi n variables: 

/>!,..., M= E fi^h,---,K) (59) 

(il ir^) 

where the sum is over all possible n-tuples (ii, ...,«„) of distinct integers 
between 1 and N. Gaudin's lemma gives a simplification of this compu- 
tation from an TV-fold integral to an n-fold integral. The Haar measure 
(see section 5.2) at the matrix A can be expressed as 

'^^=M n \e{ej)-e{ek)\^ d0i...d0N. (60) 

■ l<j<k<N 

The product here is the square of the absolute value of the N x N 
Vandermonde determinant with j, k entry e{9ky~^ = e((j — l)9k) ■ It is 
also the N x N determinant of the matrix with j, k entry JN{Oj — 9k) 
where ^ 

M9) = <rn6) = e((iV - 1)^/2)^|^ . (61) 

m=0 

Thus, 

N 

[ fiA)dA = [ m ^■■■^On)^ det Jn [Oj - e,) J] d9j . (62) 



Then, Gaudin's lemma asserts the equality 



/ f{A)dA= I f{e^^...,o^)}-AetJ^{9^-9k)\\d9 

Ju{N) J[0,1]" ~=i ■ 



(63) 



This principle works for all of the subgroups of U(i'V) under consid- 
eration here as well. (See [KS] section 5.1 for a general statement and 
proof of this important lemma.) We illustrate by computing the n-level 
density function for U(A''). Note that 

Jim lj.(./iV)=e(./2)!|^ (64) 

from which it follows easily that 

hm -J- det JNiOj - 9k) = det Kq{9i, ...,9n) (65) 

N—foo iV nxn nxn 



where is defined in section 5.5. 

Now let /(x) = f{xi, . . . ,Xn) be a test function. To compute the 
n-level density (compare with section 5.4) we need to evaluate 

lim / V f{ei„...,eijT\\e{ej)-e{ej)f de^... dON m 

By Gaudin's lemma and after using the new expression for the Haar 
measure and changing variables Oj — > Xj/N, the above is equal to 

lim / f{e,,...,0n)-,det Jn {0j -Ok) II dO, (67) 

= f{xi,...,Xn)detKo{xi,...,Xn)dx (68) 
SO that Wu,n{xi, ...,x„) = det„xn Ko{xi, a;„). 

5.5 Formulas for the Density Functions 

— The n— level density is W{xi, . . . , a;„) = det„xn K^{xi, . . . , a;„) where 
K^{xi, . . . ,Xn) is the nx n matrix with entries 

sin7r(a;i - a;^) sin7r(a;i + a;^) 
irixi-Xj) nixi+Xj) 

where e = for U; e = — 1 for Sp; e = 1 for 0+. Also, 

n 

Wo-,„(x) = det(if_i(x))+ V<5(a;„) det (ifi7^(x)) (70) 

nxn '—^ n— Ixn— 1 

m=l 

where ^ is the Dirac ^-function and the superscript m denotes that 
the m-th row and m-th column have been deleted from ii'_i(x). 

— The lowest zero density is v\ [x) where 



d 



Mx) = -^l[^^->^M) U; (71) 



^i(^)=-wrn(l- Wi(2^)) Sp; (72) 
'^i{x)=-^l[(^->^2j{2x)) O, (73) 

3=0 



where 1 > Ao(a;) > Ai(a;) . . . are the eigenvalues of 

r^/^ sin7r(<-u) 



/ rf« = A(x)/(0 



(74) 



The consecutive spacing density is 



l^ix) = Y[il~X,ix)). (75) 



— The n— correlation density, V{xi, . . . , a;„) = Wu.n{xi, ■ ■ ■ , Xn), 



5.6 The Selberg Integral 



There are many versions of Selberg's integral see |Me|; one is as follows. 
If 3?a > 0,5R/3 > 0,5R7 > -min(i f^, then 

/'•■•/' n \^^-^J\'''f[i^-^sr-H^ + ^:f-'dx, (76) 

o7n(n— 1) I I r(l + 7 + i7)^(a + J7)^(/3 + 3l) 

r(l + 7)r(a + /3 + 7(« + J-l)) ^ ^ 



5.7 Moments of Characteristic Polynomials 

Mu,n{s)= I |det(yl-/exp(-ia;))p^ (78) 
^ fr nOT.- + 2s) 



Msp,2w(s) = / I det(A - /)r (80) 

JSp(2N) 

2Ns^ + N + j)r{l/2 + s + 3 + s) 

n /^(i/2 + j-)r(i + s + iv + j) ' ^ ^ 



Mo,2n{s)= j \det{A~I)YdA (82) 

JO(2Ar) 

TV 



TT nN + J -l)r{s + J -1/2) 
f}^r{3-l/2)r{s + j+N-l)- ^ ^ 
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